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ABSTRACT 

We perform a study on the predictions of electric-dipole moments (EDMs) of neutron, 
Mercury (Hg), Thalhum (Tl), deuteron, and Radium (Ra) in the framework of next-to- 
minimal supersymmetric standard model (NMSSM) with CP-violating parameters in the 
superpotential and soft-supersymmetry-breaking sector. We confine to the case in which 
only the physical tree- level CP phase {(f)'x — <t>'^), associated with the couplings of the singlet 
terms in the superpotential and with the vacuum-expectation- values (VEVs), takes on 
a nonzero value. We found that the one-loop contributions from neutralinos are mostly 
small while the two-loop Higgs-mediated contributions of the Barr-Zee (BZ) type diagrams 
dominate. We emphasize a scenario motivated by electroweak baryogenesis. 

Keywords: Supersymmetry, Next-to-minimal supersymmetric standard model, CP violation, 
electric-dipole moments, Electroweak baryogenesis 
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1 Introduction 



Supersymmetry (SUSY) is the leading candidate for the physics beyond the standard model 
(SM). It not only solves the gauge hierarchy problem, but also provides a dynamical mech- 
anism for electroweak symmetry breaking and cosmological connections such as a natural 
candidate for the dark matter and baryogenesis. The minimal supersymmetric extension 
of the SM (MSSM) has attracted much phenomenological and theoretical interests but it 
suffers from the so-called little hierarchy problem and the fi problem. 

An extension with an extra singlet superfield, known as the next-to-minimal super- 
symmetric standard model (NMSSM) [lH6] was motivated to provide a natural solution to 
the n problem. The fi parameter in the term fiHuHd of the superpotential of the MSSM 
naturally has its value at either Mpianck or zero (due to a symmetry). However, the radia- 
tive electroweak symmetry breaking conditions require the fi parameter to be of the same 
order as the Z-boson mass for fine-tuning reasons. Such a conflict was coined as the /i 
problem [7]. In the NMSSM, the /i term is generated dynamically through the vacuum- 
expectation- value (VEV), vs, of the scalar component of the additional gauge singlet Higgs 
superfield S, which is naturally of the order of the SUSY breaking scale. Thus, an effective 
fi parameter of the order of the electroweak scale is generated. The NMSSM was recently 
revived because it was shown that it can effectively relieve the little hierarchy problem [8] . 
Due to the additional Higgs singlet field and an approximate PQ symmetry, the NMSSM 
naturally has a light pseudoscalar Higgs boson ai. It has been shown [8] that, in most 
parameter space that is natural, the SM-like Higgs boson can decay into a pair of light 
pseudoscalar bosons with a branching ratio larger than 0.7. Thus, the branching ratio of 
the SM-like Higgs boson into bb would be less than 0.3 and so the LEPII bound is effectively 
reduced to around 100 GeV [9|. Since the major decay modes of the Higgs boson are no 
longer bb, unusual search modes have been investigated [TO] . 

CP violation is one of the necessary ingredients for successful baryogenesis [11]. Al- 
though the SM can accommodate CP violation originating from the Cabibbo-Kobayashi- 
Maskawa matrix [12], it turns out that its effect is way too small to generate sufficient 
baryon asymmetry (~ 10~^°) [13] B. This fact suggests, in turn, that there should be an 
extra source of CP violation which has not been probed yet. CP violation relevant to 
electroweak baryogenesis (EWBG) [15] by construction must appear in the Higgs self in- 
teractions and/or the Higgs interactions with the other particles whose masses are 0(100) 
GeV. Therefore, such CP violating effects can be communicated to the low energy observ- 
ables which are measurable in the near future experiments. A lot of effort on the EWBG 
study have been made in the new physics models such as the MSSM [16], the two-Higgs 
doublet model [17] and the singlet-extended MSSM [T8l[T9]. 

* Another shortcoming in the SM baryogenesis is that the electroweak phase transition is a smooth 
crossover for mh > 73 GeV [14;, rendering thermal nonequilibrium unrealizable. 
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The MSSM offers many possible sources of CP violation beyond the single Kobayashi- 
Maskawa phase in the SM. As far as the Higgs sector is concerned, the non-vanishing CP 
phases could induce significant mixing between the CP-even and CP-odd states radia- 
tively [201423] , giving rise to a number of interesting CP violating phenomena and substan- 
tial modifications to Higgs-boson phenomenology [2U[25]. In particular, the lightest Higgs 
boson can be as light as a few GeV with almost vanishing couplings to the weak gauge 
bosons when the CP- violating phases are maximal. The decay patterns of the heavier Higgs 
bosons become much more complicated compared to the CP-conserving case because of the 
loss of its CP parities [26l|27]. These combined features make the Higgs boson searches at 
LEP difficult, consequently, the Higgs boson lighter than ~ 50 GeV can survive the LEP 
hmit [2H]. 

The non-observation of electric dipole moments (EDMs) for Thallium [29], neu- 
tron [3U], and Mercury [211122] is known to constrain the CP- violating phases very tightly. 
It is generally believed that one- loop contributions dominate and we set the phases to O{0) 
to make the null results of the EDM searches consistent within most of the parameter space 
However, we point out in this work that it may not be the case in the framework of 
NMSSM with CP-violating parameters. Even if we set the CP phases of the parameters 
appearing in the MSSM to zero, there could be potentially large nontrivial two-loop contri- 
butions coming from a combination of the CP phases, {(j)'^ — (p'lJ, which could exist only in 
the NMSSM. We note that, being different from the MSSM, the non-vanishing CP phase 
could cause CP-violating mixing among the neutral Higgs bosons even at the tree level. 

In this work, we perform a study on the predictions of EDMs of neutron. Mercury 
(Hg), Thallium (Tl), deuteron, and Radium (Ra) in the framework of NMSSM with CP- 
violation. We confine ourselves to the case in which only the physical CP phase {(j)')^ — 0'^), 
associated with the couplings of the terms containing the singlet field in the superpotential 
and with the VEVs, takes on a nonzero value. We figure out how large the CP phase 
can be taken in a scenario in which a first-order phase transition could be achieved more 
easily in comparison to the MSSM [IB]. The form factors that contribute to these ob- 
servable EDMs include electric-dipole moment, chromo-electric dipole moment, Weinberg 
three-gluon operator, and the four-fermion operators. The two-loop Weinberg three-gluon 
operator and the Higgs-exchange four-fermion operators are generated due to the tree- 
level CP-violating Higgs mixing. The electric-dipole moment (EDM) and chromo-electric- 
dipole moment (CEDM) receive the following one- and two-loop contributions: (i) One 
loop neutralino-sfermion contribution in which the CP phase appears in the neutralino 
mass matrix (the CP phase of the effective /i parameter in the chargino-mass matrix is set 
to zero), (ii) Two- loop Barr-Zee (BZ) diagram with the 7iJ°, W^H^, W^W^, and ZH^ 
decompositions. 

^Nevertheless, accidental cancellations among various contributions may occur in the three measured 
EDMs, thus still allowing sizable CP phases even with the SUSY particles lighter than 0(1 TeV) 
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We found that the one-loop contributions from the neutrahno-sfermion diagrams to 
the Thalhum and neutron EDMs he below the present experimental upper limits especially 
when the sfermions of the first two generations are heavier than ~ 300 GeV. This is in 
agreement with the previous observations |11[35], in which only the one- loop contributions 
were taken into account. The one-loop contributions to the Mercury EDM could be larger 
but they also go below the present experimental upper limit if the sfermions of the first 
two generations are heavier than ~ 300 — 500 GeV. 

The two-loop contributions start to dominate when the sfermions of the first two 
generations are heavier than ~ 300 GeV and the one-loop contributions are suppressed. We 
found that the two-loop contributions can saturate the current bound on the neutron EDM 
and they can go over that on the Mercury EDM. But we found that there is still a room 
to have the maximal CP phase {cj)')^ — 0'^) ~ 90° after taking account of the uncertainties 
in the calculations of the EDMs. We note that the large CP phase can be easily probed 
in the proposed future experiments searching for the EDMs of the deuteron and the ^^^Ra 
atom and it might be connected to the EWBG. 

The organization of the paper is as follows. We briefly describe the Higgs sector of 
the NMSSM with CP-violating parameters in Sec. II. We give the relevant Higgs couplings 
in Sec. Ill and detail breakdowns of the EDM calculations in Sec. IV. Numerical analysis 
is given in Sec. V. We conclude in Sec. VI. 



2 Higgs sector in the NMSSM with CP violation 

The superpotential of the NMSSM may be written as 

H^NMSSM = U'^KQEu + D^'YidHaQ + E^'KHoL + XSHJId + (1) 

where S denotes the singlet Higgs superfield, H^^ are the two SU(2)i doublet Higgs super- 
fields, and Q, L and , D'-^, E'^ are the matter doublet and singlet superfields, respectively, 
related to up- and down- type quarks and charged leptons. We note that, especially, the 
last cubic term with a dimensionless coupling n respects an extra discrete Z3 symmetry. 
The superpotential leads to the tree-level Higgs potential, which is given by the sum 

Vo = Vf + Vd + Koft, (2) 

where each term is given by 

Vf = \\\^\S\\HlHd + HlHu) + \\HuHa + KS^\\ 

Koft = mlHlH, + mlHiH^ + ml\S\^+(xA^SH^H,-l^A,S^ + h.c), (3) 



4 



with the gauge-couphng constants g' = e/ cos6'vy and g = e/ sin6'vy, where e is the electric 
charge of the positron in our convention. Note that we are taking the unusual minus(— ) 
sign for the singlet soft-trilinear term proportional to A^. 

We have parametrized the component fields of the two doublet and one singlet scalar 
Higgs fields and the vacuum expectation values (VEVs) as follows, 



S = —{vs + <j^"s + ^as). (4) 

Note that we have the complex vacuum-expectation-values (VEVs) and assume that the 
parameters A and k in the superpotential and Ax and in the soft terms contain non- 
trivial CP phases. It turns out that not all the CP phases appearing at the tree level after 
the electroweak symmetry breaking are physical and the only physical one is the difference 
- <P'k with 

</'a = </'A + 6* + V? and (p'^ = + Syp , (5) 

and the CP phases of Ax and A^ are determined up to a two-fold ambiguity using the two 
CP-odd tadpole conditions. When (p'x — 0'^ 7^ 0, the neutral Higgs bosons do not have to 
carry any definite CP parities already at the tree level and its mixing is described by the 
orthogonal 5x5 matrix O^i as 

°, 00,0°, a, 05)^ = 0^i{Hi,H2,Hs,H,,H,f (6) 

with -ffi{5) the lightest (heaviest) Higgs mass eigenstate. 

For the calculation of the Higgs-boson masses and mixing matrix Oai in the presence 
of CP-violating parameters in the superpotential and in the soft-supersymmetry-breaking 
sector, we adopt the renormalization-group (RG) improved approach by including the full 
one-loop and the logarithmically enhanced two-loop effects |36]. And then, the NMSSM 
Higgs sector is fixed by specifying the following input parameters: 



tree level 
1-loop level 
CP phases 
signs of 



|A| , |k| ,tan/3 ; \Ax\ , \A^\ ,vs 
Mq.^,M^^,Mj^^,\A,\,\A,\ 

<P'\ , ; , 0^6 

Lf < J. \ / aI 



COs{(j)'x + 0A J , COS(0'^ + 0A J • (7) 

For the renormalization scale Qq we take the top-quark mass as in Refs. 
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In this work, we wish to consider the constraint on the tree-level CP phase com- 
ing from the non-observation of electric dipole moments (EDMs) for Thallium (^°^T1) [29] . 
the neutron (n) [30], and Mercury (^^^Hg) [3T1|32] when the CP phases appearing in all the 
other soft SUSY-breaking terms vanish or sin(</)';^ + (pAf) = sin(0';^ + (pi) = 0, with (p^^. and 
(pi denoting the CP phases of the soft trilinear parameters Af and the three gaugino mass 
parameters Mj=i^2,3, respectively. 



3 Higgs-boson couplings in the NMSSM 

If we consider the case in which only the tree-level CP phase (p'^ — (p'^ takes a non-trivial 
value while all the other CP phases are vanishing, as will be shown in the following, the 
one-loop contributions to the EDMs are mostly small and the EDMs are dominated by 
the two-loop contributions from the Higgs-mediated dimension-6 Weinberg operator, the 
Higgs-exchange four-fermion operators, and the Barr-Zee-type diagrams. Therefore, for 
the calculation of the EDMs beyond the one-loop level, one may need the Higgs-boson 
couplings taking full account of the 5x5 CP- violating mixing matrix Oai- In this section, we 
present the couplings of the neutral and charged Higgs bosons to quarks, leptons, charginos, 
neutralinos, and third-generation sfermions in the NMSSM with CP violation. For the 
conventions and notations of the masses and mixing matrices of charginos, neutralinos, 
and third-generation sfermions, we refer to Appendix A. 

The interactions of the five neutral Higgs bosons with the SM quarks and leptons are 
described by the interaction Lagrangian: 
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i=l 



where gf = gmf/2Mw for f = u,d, I. At the tree level, {gjj.ff, gn ff) ~ (^i*/"^/?' ~Oii tan/3) 
and {guff' duff) ~ {02i/ sp, —Oacot (3) for / = (/,(i) and f = u, respectively. The 
simultaneous existence of the couplings gfj.ff and g^ jf signals CP violation. The couplings 
of the neutral Higgs bosons to the five neutralinos are given by: 



^i/o~o~o = -fE^fcX? (4,^0-0 +275<5;ro~o)x; (9) 

where i, j = 1, 2, 3, 4, 5 for the five neutralinos and A; = 1, 2, 3, 4, 5 for the five neutral Higgs 
bosons and the scalar and pseudo-scalar coupling are 

akxlx^ = 3f?e|^(A72-%Ar;,)[Ar^,(Oi,-204fcS^)-A^;^(02fc-204fcC^)] 

(Osfc + iO,,)n:,n;^ + {02k + iO,kCp)N:,N;^ 



V2 



g 
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-{Oik + iOikSp)N:^N. 



J4 



|^|g«(</'K+3ip) 



V2 
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{0,u + iO,u)N*,N*, +[i^3] 



|^|g«(</'K+3v') 



V2 
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{0^k + iO,u)N:,N*\+[i^j]\. 



(10) 



Note the couplings are symmetric under the exchange oi i O j, reflecting the Majorana 
property of the neutralinos, and contain the terms coupled to the singlet components of 
the Higgs bosons and to the singlino components of the neutralinos which do not exist in 
the MSSM. The couplings of the neutral Higgs bosons to the charginos can be similarly 
cast into the form: 



'ID 



with i,j = 1, 2 for the two charginos, = 1, 2, 3, 4, 5 for the five neutral Higgses, and 



^HkxtXj 
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HkxtXj 



+ (Ci?,)*2(CL)*2(03fc + ^Osfc) + [l ^ ]]* \ , 

{CRh{CL)*2{0ik - ^0,kSfs) + {CR),2{CL)*i{02k - lO^kCp] 
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-{CR),2{CL)UO-ik + tO,k)-[t^3r). 



(12) 



We observe the couplings are real when i = j but they are complex otherwise with 



9h, 



S,P 

kxtxT 



S,P 
kxtx2 



3h 



We again note that the couplings contain the terms coupled to the 
singlet components of the neutral Higgs bosons which do not exist in the MSSM. Lastly, 
the neutral Higgs interactions with sfermions can be written in terms of the sfermion mass 
eigenstates as 

(13) 
(14) 



Hff _ ... 

f=u,d,l,u 



where 



vg 
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with a = (00, 4, a, as) = {l, 2, 3, 4, 5) , /3 , 7 = (/l, /«) = (1, 2) , z = (ifi, ifa, i^s, H^,H^) 
(1, 2, 3, 4, 5) , and j ,k = (/i, = (1, 2). The exphcit expressions of the couphngs r^^'^ in 
the weak basis are given in Appendix B. 

Now let us move to the couphngs of the charged Higgs bosons. The charged Higgs 
boson couphngs to the SM quarks and leptons are described by the Lagrangian 

(15) 



with gf^f^ = —grriu/ ^f^mw and —gmi/ yj^mw when {f^, fi) = {u,d) and (uj), respec- 
tively. At the tree level, {gH+ud ^ 9H+ud) = {[l/tp + {md/mu)tp)]/2 ,i[l/tp- {md/mu)tp)]/2) 
and idH+pi y 9h+pi) = p / ^ p / ■ The interactions of the charged Higgs bosons with 
charginos and neutralinos are described by the following Lagrangian: 

9 



with 2 = 1,2, 3, 4, 5, j = 1, 2, and 



(16) 



X.i Xj 



i2 



C/3 



+ 



V2\X\ 



R}j2 



p 

9h+y0 



Xi Xj 



V2\\\ 



L)j2 



C/3 



V2\X\t 
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(17) 



As similarly in the neutral Higgs couplings, we note that the couplings are containing, 
in addition to the corresponding MSSM interactions, the terms coupled to the singlino 
components of the neutralinos. Finally, the charged Higgs couplings to sfermions are given 
by 



H±f'f 



vg 



H+flf'k 



H+fU\] +h.c. 



where 



vg 



(19) 



The explicit expressions of the couplings T^"^^*^' in the weak basis are also given in Ap- 
pendix B. 
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4 Synopsis of EDMs 



In this section, we briefly outline how we estimate observable EDMs. We start by giving 
the relevant interaction Lagrangian as follows: 



L = - ^ rff F^^ / a^,75 / - ^ < G"^ '''^ q a,,^,T^ q 



I 

+\d'' U,G';,^&^''G\^ + EC//'(/7)(/W) , (20) 

where F^'^ and are the electromagnetic and strong field strengths, respectively, the 

rj^a ^ ^^12 are the generators of the SU(3)c group and G^'' = lei'^'^'^Gxa is the dual of the 
SU(3)c field-strength tensor Gao-- The EDMs and the CEDMs of quarks and leptons are 
denoted by and d^, respectively. 

For the Weinberg operator, we consider the contributions from the Higgs-mediated 
two-loop diagrams: 



Wj ,=t,b 



dUiqq 9H,qq ^{Ziq) 



(21) 



where Ziq = Mfjjrn?^ and, for the loop function h{ziq), we refer to Ref. [38] . 

For the four-fermion operators, we consider the t-channel exchanges of the CP- 
violating neutral Higgs bosons which give rise to the CP-odd coefficients as follows [M] : 

(Cirf = SfSf'j: '"-"'"'"' ■ (22) 

i Hi 

The EDM dj and CEDM d^' are given by the sum of the one-loop and two-loop 
contributions 

df = {dff + {dfr ; < = «)^° + (d^r ■ (23) 

The details of the neutralino-mediated one-loop contributions and the contributions from 
the two-loop Barr-Zee-type diagrams will be discussed below. 



4.1 One-loop EDMs 

In the case under consideration, the only non-vanishing one-loop contribution to the (C)EDMs 
comes from the neutralino loops due to the CP phase (f)'^ — 0'^. The one-loop contributions 
to the EDMs of charged leptons {df /e)^^\ up-type quarks {df/e)^° and down- type quarks 
{df /e)^" may conveniently be expressed as 

1 5 2 — — — 

E E -1^ ^HiaiY/r 9{ff] QjB{mym)) , (24) 




IGvr^ ^ m~ 
1=1 j=i J, 
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with f = l,u, d. The neutrahno-fermion-sfermion couphngs are 



ij — V^y^3iVj2V'^ ;ij - V I'M/ VV/ -^3 ^^^ilV^ Jlj ~ "'I^'iaK'^ )2j ^ 

g^i/ = V2gtwQfNa{U%-h}NUU%, (25) 



where the Higgsino index a = 3 {f = l,d) or 4 (/ = u), T^'*^ = — 1/2 and = +1/2 and 
the loop function is given by 

with -B(l) = 1/6. As well as the EDMs, the neutralino loops can induce non- vanishing 
chromo-electric dipole moments (CEDMs) for the quarks as follows: 

{dU^f = if^ EE 5^ ^-[(,i:f)*,?,-]5(rn|o/rn|). (27) 
-L"^ i=ij=i 

4.2 Two-loop Barr-Zee EDMs 

Beyond the one-loop, we take account of the contributions from the two-loop Barr-Zee- 
type diagrams. We have considered the the Barr-Zee diagrams mediated by the ^-'y-H^ 
couplings [31] and the ■y-H'^-W^ and 'j-W'^-W^ couplings [39]. The two-loop diagrams 
mediated by the '-/-H^-Z couplings [40ll4T] have also been included taking account of the 
general CP-violating Higgs-boson mixing. More explicitly, the contribution from the two- 
loop Higgs-mediated Barr-Zee-type diagrams can be decomposed into four parts: 

(df) = {dfY + {df) + (df) + (df) (28) 

where 

, O r- X ('''Y"' - V f^2;;iOS!!Vf <,7 V „ Fir ) 



I ^Q^cm Qq ^/ \ P _ S f( \ \ S _ P ( ^ 
Svr^S^M^ ^ [^H,/75'/f,gq / l^gi J + 9Hjf9H,gq9{'rqi, 

aemrrif gjjj^ V o ~ ~ Fir- 1 

1=1 rli J = l,2 

^ H-^^T^ M2 E [9Hjf9H,r+r- f i^ri) + gHjf9H,r+r- 9{rri] 
"em 
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5 2 

xEE 



a 



w 



1 



1 — X 

s 



J rwH± 



Ml, 



r 



: ^X'lH±' 



dx 

+ ^g^^,J^f;) G^J") m~o (1 - x) 



1 — X X 
P \ r^RL^ 



'd 



a 



5 2 



EE 



327r2s^ i=i j=i 



^VF+x°x7 I V+x?> 



(29) 



(30) 
(31) 



with r,, = ml/Ml^, r^y = M^/M^ rj = m^/M^ and r, = m^jMl, and (/t, /J 

A. J A j 

(n, (i) , (z/; , /). The l^-boson couphngs to the charginos and neutrahnos are given by 



iV.3(CL)*2 + V^^^^2(CL)*l , 



and, with A,B = L, R, 



w+x"x" 



'w+x1x~ 



(32) 



.(33) 



For the loop functions -F(r), /(t), g{r), J{a,b), and /vi/VK('"i , '"j), we refer to, for example, 
Refs. [M|[39] and references therein. Finally, for (dj )^^°, we take account of the dominant 
fermionic contributions given by 



'd 



a V 



zff 



E 



SQgTUq 



X 



E 

i=l 



rrir. 



m% Jo 



1 1 / 
dx-J \ rzH, 
X \ 



~3Hjf ( 



VZqqQUiqq 



m% JO 



dx J rzHi, 



x{l — x) J 

rqH, 



X 



x{l — x) ^ 
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a Vzff rrif rrij 
16V27r2cLsL ^ 72M, 



w 



a'^Vzff TUf 



mjj, Jo X 



1 1 / 
dx-J rzH,, 



iw 



5 2 

X EE 

k=li,j=l L 



mi. JO 



5 



dx- -J [rzH, 



x{l — x) ) 



X 



x{l — x) ^ 



X 



i ^kXj^ Xj 



1 1 
dx-J I rzH, 



^^Zx+x79H^xtx~ 



xr~~rr + (1 — x)r 

Xi k 



x7Hk 



k ' 



x(l — x) 



+9kff^'^('''zxtx-9H, 



_ _ p 
^Zx+xT^HkxtxJ 



m~ 



X 



1 1 -X 
dx 



xr~-fj + (1 - x)r; 



XJH, 



x{l — x) 



X, '^zxjx- 



75 ) Xj Z^, 



'Zx7x- 



The Z-boson couplings to the charginos are given by 

^zx+x- = -azxTi" [vzxt-y 

where gz = g/cw = e/{swCw) and 

And the Z-boson couphngs to the quarks and leptons are given by 
with Vzff = Tij2-Qfsl, and 

<^zff = tIl/2. For the SM quarks and leptons, T^^'' 



'Zxtx- 



(34) 



(35) 



(36) 

(37) 
+ 1/2 

In addition to EDMs, the two-loop Higgs-mediated Barr-Zee graphs also generate 
CEDMs of the light quarks qi = u,d which take the forms: 



and T-: 



d,e 



BZ 



-E' 

q=t,b 



647r3 Ml . 

1=1 tli j = l,2 



+ 



gs ^em. 



E [9H,qiqi9Hiqq f i'^qi) + 9H,qiq,9H,qq 9{'^qi] 



{31 
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4.3 Observable EDMs 



In this subsection, we briefly review the dependence of the Thallium, neutron. Mercury, 
deuteron, and Radium EDMs on the (C)EDMs of quarks and leptons and the coefficients 
of the dimension-six Weinberg operator and the four-fermion operators. 



4.3.1 Thallium EDM 



The Thallium EDM receives contributions mainly from two terms [l2l 

dTi [e cm] = -585 ■ rff [e cm] - 8.5 x 10"^^ [e cm] ■ (Cs TeV^) + ■ ■ • , (39) 

where is the electron EDM and is the coefficient of the CP-odd electro n-nucleon 
interaction Ccs = e^75 eNN which is given by 

, ^C. r + (0.1 GeV) ^ t %^ (40) 

with n = {N\m,ss\N) /220 MeV ~ 0.50 ± 0.25 and 



q=t,b [ j=l,2 "^qj 



with Xt = 1 and Xb = 1 — 0.25/t. 



4.3.2 Neutron EDM 



For the neutron EDM, we consider three different hadronic approaches: (i) the Chiral 
Quark Model (CQM), (ii) the Parton Quark Model (PQM) and (iii) the QCD sum-rule 
technique. 



In the CQM approach, the neutron EDM is given by 

^ _ l^NDA _ 1 jNDA 

3 3 ' 

Cm = ^^< + ^^i< + ^^^^^ (42) 

where the chiral symmetry breaking scale A ~ 1.19 GeV and the rj^''-''^ account for 
the renormalization-group (RG) evolution of dg''" and d^ from the electroweak (EW) 
scale to the hadronic scale. For the QCD correction factors we are taking r]^ ~ 1.53 
and 77*^ ~ r/*^ ^ 3.4 [H]. We note that the EDM operators d^''^ and d*^ in f H2]) are 
computed at the EW scale. 
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In the PQM approach, the neutron EDM is given by [15] 



(43) 



with A^«^^ 



0.746, APQ^ = -0.508, and AfQ^ = -0.226. 
Using the QCD sum-rule technique, the neutron EDM is given by [lEHSO] 



dn{d^ ,d'^) + dnid"^) + dniCbd) + 



dn{d 



G\ 



dn{C) 



bdj 



(1.4 ±0.6) (rff -0.25df 
±e(20± 10) MeVd^, 

±e2.6 X 10"^ GeV^ 



;i.l±0.5)e(d^ + 0.5rf:;)/^?. 



+ 0.75 

7716. 



(44) 



where c?^ and d'^ should be evaluated at the EW scale and d"-^ at the 1 GeV scale, 

Bl. In the numerical estimates we take the 



d<^\. 



GeV 



(r/^/0.4) d^l 



EW 



\.bd^\ 



EW 



positive signs for dn{d ) and dn{Cbd)- 



4.3.3 Mercury EDM 



Using the QCD sum rules |l9l|50], we estimate the Mercury EDM as follows: 



d 



I, II ,111, IV 



[S] + 10^'C + (3-5 X lO-'GeV) eCs 
Z -N' 



+ (4 X 10-* GeV) e 



Cp + 



A 



C' 



(45) 



where d 



1 ,11 ,111 ,iv 

Hg 



[5*] denotes the Mercury EDM induced by the Schiff moment. The pa- 



rameters Cp and C'p are the couplings of electro n-nucleon interactions as in Ccp 
Cp ee Ni'j^N + C'peeNi'j^TsN and they are given by 



Cr 



Op _ 



a. 



375 MeV J2 — 



q=c,s,t,b y 



806 MeV 



ed 



md 



- 181 MeV 



C, 



eq 



q=c,s,t,b y 



(46) 



We take account of the uncertainties in the calculation of the Schiff-moment induced Mer- 
cury EDM as follows [51] : 



d^lS] 



1.8x 10-=^e^S^/GeV, 

7.6 X 10-^e^S^/GeV + l.O x lO'^ e^S^ /GeV , 
1.3x 10"*e^S^/GeV + 1.4x lO-^e^S^/GeV, 
3.1 X 10"* e gflj^ /GeV + 9.5 x 10"^ e /GeV . 



(47) 
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where 



-(0) 
9nNN 



y-nNN 



0.4 X 10 



2t\ X 10 



-12 



{d'^+(f^)/gs m)\ 



10-26cm (225MeV)3 

-12 (rf? - rf.^)/^?. m)\ 



10-26cm (225MeV)3 
8 X 10-=^GeV=^ \^-^^ + 3.3/€^ 4 



'I - 0.25ft;) 



bd 

rrib. 



(4^ 



4.3.4 Deuteron EDM 



For the deuteron EDM, we use [Ml 152] : 



d 



D 



-11 



(0.6 ±0.3)1 e(rf^-rf^)/^?. 

+ (( 

0.5a 



-(0.2±0.1)e(d^ + rf^)/<7, + (0.5±0.3)Kf + rff) 



+ (1 ±0.2) X 10^2 eGeV^ 
± e(20± 10) MeVd^ . 



3.3^^ + (1 



0.25k) 



C, 



bd 



nib 



(49) 



In the above, d'^ is evaluated at the 1 GeV scale, and the coupling coefficients gd,s,b ap- 
pearing in Cdd,sd,bd are computed at energies 1 GeV, 1 GeV and nib, respectively. All other 
EDM operators are calculated at the EW scale. In the numerical estimates we take the 
positive sign for d*-^. 



4.3.5 Radium EDM 

For the EDM of ^^SRa, we use [51]: 

rfRa ^ dR.4S] - -8.7 X 10-2 e^S^ /GeV ± 3.5 x 10"^ e /GeV . (50) 

We note that the ^^^^^v contribution to the Radium EDM is about 200 times larger than 
that to the Mercury EDM dL[S] [53]. 



5 Numerical Analysis 

The scenario we are considering has an intermediate value of tan/3 with small vs ~ v: 
tan/3 = 5, ^75 = 200 GeV, 

Mp^ =Mrj =M^ =Mj =M^ =1 TeV , 

VI, 2, 3 f^l,2,3 ^1,2,3 -'^1,2,3 ^1,2,3 

l^el = \A^\ = \Ad\ = \As\ = \A^\ = \At\ = \Ab\ = 1 TeV, 

0Ae = (pA^ = (pAa = <PAs = <pAr = <pAt = = ; 0l,2 = TT , 03 = , 

= ; sign [cos(0'^ ± 0aJ] = sign [cos(0';, ± 0^J] = ±1 , (51) 
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Figure 1: The absolute value of the Thallium EDM d^i divided by the current experimental 
limit d^f^ as a function of 0'^ varying |A|, \Ax\, and l^^l over the ranges given by 
Eq. (l55ll for the scenario specified by Eq. ( ISTl) . Especially, the |k| range is divided into 2 
regions: 0.01 < |k| < 0.1 (red filled circle), and 0.1 < |k| < 0.2 (blue plus) . 

while varying 

|A|,|k|; \Ax\,\A^\; (j)'^. (52) 

We have chosen Mi = M2 = —200 GeV to fix the neutralino sector 0. We find that 
a first-order phase transition could occur in some regions of the parameter space of this 

^In some regions of the parameter space, we find that T{Z — >• XiXi) becomes sizable around (j)'^ = 180°, 
violating the LEP bound on the non-SM contributions to the invisible Z decay width, AFinv < 2 MeV [51] . 
In this section, we are presenting our results without including the bound on AFinv since it can be easily 
satisfied for other choices of Mi^2 without affecting the numerical results much. 
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Figure 2: The same as in Fig. [T]but for the neutron EDM in the QCD sum rules (left), 
CQM (middle) and PQM approaches (right). 

scenario [18] which is needed for the EWBG [15] . 

In Figs, m El and El we show the absolute values of 

d^i/ c^ff ^ 1 , c?Hg/ c^Hg^ ) (53) 

with the current experimental bounds 

cilfP = 9 X 10-25ecm, rf^^^ = 2.9 x lO"'^ e cm , = 3.1 x lO"'^ e cm . (54) 

For a given value of 0'^, we perform a scan by sampling the four model parameters in the 
following ranges: 

0.75 < |A| <0.95, 
0.01 < \k\ <0.2, 
|Aa|min< \Ax\ < 800 GeV, 

|A|min< I^.I < 200 GeV, (55) 

where |Aa,k|min are determined by the tadpole conditions [36]. The lower limit on |A| 
comes from the chargino mass limit and that on \k\ is derived from the global minimum 
condition and the requirement of the strong enough first-order electroweak phase transition. 
The upper limits on |A| and \k\ come from requiring that there is no serious breakdown of 
perturbativity below the GUT scale. Especially, we have taken 0.95 as the upper bound for 
|A| as done in Ref. [18], though it is somewhat larger than the usual perturbativity bound 
(~ 0.8) quoted in the NMSSM. Such a value of |A| might be comfortably accommodated 
in, for example, the SUSY fat Higgs model which includes a new gauge interaction that 
becomes strong at an intermediate scale [55] . And then the LEP limits, the global minimum 
condition, and the positivity of the square of the Higgs-boson mass have been imposed as 
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Figure 3: The same as in Fig. [T]but for the Mercury EDM: d^^ (upper left), (i^g (upper 
right), rf^g (lower left), and rf^g (lower right). 

described in Ref. [36]. We note that the region of \Ax\ is around |A|f5tan/3/A/2 ~ 600 
GeV, which determines the typical masses of the heavier Higgs bosons. 

The Thallium EDM is below the current experimental limits over the whole range of 
the parameters except for a few points around 0'^ = 110° and 250°, see Fig. [1] Also, we see 
that the ratio \dTi/d^f-^\ does not exceed 3. 

Fig. |2] shows the neutron EDM in the three different approaches. The estimations in 
the QCD sum rules and CQM approaches give more or less similar results which do not 
exceed 3 times the current experimental limit while the PQM prediction always lies below 
it. We observe that larger values of \k,\ lead to larger EDMs in the QCD and PQM cases. 
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Figure 4: The same as in Fig. [T] but for the deuteron (left) and Radium (right) EDMs. 

Fig. [3] shows the four different calculations of the Mercury EDM. Using the three 
calculations d^^, d^^, and d^^, we observe that I^Hg/c^Hg^l can be as large as ~ 20 and so 
only the small angle region with A0'^ ~ ±10° (30°) around 0° is allowed when |k| > (<)0.1. 
The small angle region around 180° is also allowed. While if we adopt d^^, the Mercury 
EDM does not exceed 4 times the current experimental limit and much larger values of 0'^ 
are still allowed. 

To summarize, the tree-level CP phase 0'^ is hardly constrained by the non-observation 
of Thallium and neutron EDMs. The Mercury EDM constraint could be stronger but there 
is still a room to have large (/)'^ ~ 90° after taking account of the uncertainties in the 
calculation of the Schiff- moment induced Mercury EDM, dngi'S']. 

At this stage, it would be interesting to see whether the proposed future EDM ex- 
periments can probe the CP phase (/)'^. In this work, we consider the deuteron EDM and 
the ^^^Ra EDM. The latter is known to be enhanced by the Schiff moment induced by the 
presence of nearby parity-doublet states [53]. In Fig. HI we show the absolute values of 

d^/dl'''' and dnJdlY. (56) 

For the normalization of the deuteron EDM, we have taken the projected experimental 
sensitivity [56]: d^^ = 3 x 10~^^ecm. For the Radium EDM, we have taken d^^ = 
1 X 10~^^ecm. The chosen value for d^^ is near to a sensitivity which can be achieved in 
one day of data-taking [57]. From Fig. HI we see that I^d/'^d'^^I MRa/c^Ra^l can be as 
large as ~ 150 and the larger values of \k,\ lead to the larger EDMs. We observe that both 
of the experiments have powerful potential to probe almost the whole region of 0'^. 
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50 100 150 200 250 300 350 

([);^ [degrees] 

Figure 5: The observable EDMs taking |A| = 0.81, \k\ = 0.08, \Ax\ = 575 GeV, and 
\Ai^\ = 110 GeV. The other parameters are fixed as in Eq. (15T1) 

In the remaining part of this section, we wish to present the details of the observable 
EDMs by exemplifying a point 

|A| = 0.81, |k|=0.08, \Ax\ = 575 GeV, = 110 GeV. (57) 

The other parameters are fixed as in Eq. (|5T1) . as motivated by EWBG. In Fig. [5l we show 
the absolute values of the observable EDMs under consideration divided by the correspond- 
ing current experimental limits or the projected experimental sensitivities. 

Fig. [6] shows the Thallium, neutron and Mercury EDMs together with the constituent 
contributions. We observe that both the df and Cs terms significantly contribute to the 
Thallium EDM. The dominant contributions to the neutron EDM in the QCD sum-rule 
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Figure 6: The observable EDMs together with the constituent contributions taking |A| = 
0.81, |k| = 0.08, I^aI = 575 GeV, and \Ai^\ = 110 GeV. The other parameters are fixed as 
in Eq. ( ISTl) . Each frame shows: the Thalhum EDM (upper left), the neutron EDM in the 
QCD sum- rule approach (upper right), the neutron EDM in the PQM (lower left), and the 
Mercury EDM d^^ (lower right). 

approach come from the CEDM and d'^ terms and we note that the neutron EDM in the 
CQM shows the similar behavior (not shown). The neutron EDM in the PQM is dominated 
by the contributions from the EDMs of the up and strange quarks. On the other hand, the 
Mercury EDM is dominated by the CEDMs of the light quarks. Fig. [7] shows the deuteron 
and Radium EDMs together with the constituent contributions. We observe both of the 
EDMs are dominated by the contributions from the CEDM terms. 
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Figure 7: The observable EDMs together with the constituent contributions taking |A| = 
0.81, = 0.08, \Ax\ = 575 GeV, and \A^\ = 110 GeV. The other parameters are fixed as 
in Eq. flSTj) . The left (right) frame shows the deuteron (Radium) EDM. 

Fig. |8] shows the EDMs and CEDMS of the electron and light quarks together with 
the constituent contributions. We observe that the electron EDM is dominated by the two- 
loop Barr-Zee contributions mediated by the 7-7- if? and j-W^-W^ couplings, whereas the 
one-loop contribution from the neutralino loops is subleading. The one-loop contribution 
is suppressed because the CP phase (j)'^ can contribute to EDM only through the multiple 
singlino-Higgsino-gaugino mixing and we are taking somewhat large values for the masses 
of the sfermions of the first two generations: Mp^ = M^, = = Mr = M9; = 1 

Ql,2 '-'1,2 Ul,2 1^1,2 1^1,2 

TeV. The other two-loop contributions mediated by the '^-H^-W^ and '-/-H^-Z couplings 
are suppressed by the large charged Higgs-boson mass ~ 600 GeV and the small vector 
coupling of the Z boson to electrons, Vze+e- = + with = 0.23, respectively. 

The EDMs of the light quarks are dominated by the two-loop Barr-Zee contributions 
mediated by the 'j-j-Hf and 'j-W^-W^ couplings. Being different from the case of d^, 
the two-loop contribution mediated by the 'y-H^-Z couplings is larger than the one-loop 
contribution. We note a cancellation occurs between the two dominant two-loop Barr-Zee 
contributions in d^^, resulting in \d^\ < \d^\ ~ \df\. The CEDMs of the up and down 
quarks are dominated by the Higgs-mediated two-loop contributions which are more than 
100 times larger than the one-loop contributions and we note l/M^ | ~ {rnd/rriu) tan/3 ~ 
10. 

Finally, we examine the dependence of the observable EDMs on the sfermion masses 
of the first two generations. Fig. [9] shows the observables EDMs as functions of the common 
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mass scale for the first two venerations, = Mp; = M^r = M-x = Mr = . 

' ^1,2 Ql,2 Ul,2 Ul,2 1^1,2 -C/1,2 

Except the neutron EDM based on the CQM which hes below the current experimental 
limit independently of ^ , all the EDMs exhibit dips at certain values of ^ . The 
dips occur because of the cancellation between the one- and two-loop contributions. When 
M^^ ^ is small the one-loop contribution dominates. As M^^ ^ grows, the one-loop con- 
tribution decouples and the EDMs saturate to certain values determined by the two-loop 
contribution. Therefore, for the neutron (QCD), Mercury, deuteron, and Radium EDMs, 
we observe that the one-loop contribution is comparable to or larger than the two-loop one 
only when ^ < 300 GeV. On the other hand, for the neutron EDM based on the PQM 
and the Thallium EDM, the one-loop contributions are larger but the two-loop contribution 
starts to dominate when M^^ ^ is larger than ~ 600 GeV. By choosing for the Mercury 
EDM, we observe that all the EDM constraints could be fulfilled when M~ > 300 GeV 
without relying on the cancellation mechanism. If we make other choices for the Mercury 
EDM, to suppress all the EDMs for Thallium, neutron, and Mercury below their present 
experimental bounds, the required degree of cancellation is about 90 % over the whole 
range of with 100 % corresponding to complete cancellation. 

Before we close this section, we make a comment on the tan /3 dependence of the EDMs 
of the electron and the down and strange quarks. The one-loop neutralino contribution is 
proportional to tan/3. The Barr-Zee contribution mediated by the j-W^-W^ couplings, 
one of the two two-loop leading contributions, is nearly independent of tan/3. The tan/3 
dependence of the other dominant Barr-Zee contribution mediated by the 7-7- if? couplings 
is much milder compared to the MSSM case in which the sbottom and stau contributions 
are proportional to tan^/3 [IHlESllSn]- This is because the masses of the two heavy Higgs 
states, which include the CP-odd state from the Higgs doublets, increase as tan /3 grows, 
Mh^^Hs — ilM'^^s/V^) tan/3, to avoid tachyonic Higgs states [36] . 



6 Conclusions 

We have performed a study on the predictions of EDMs for Thallium, neutron. Mercury, 
deuteron, and Radium in the framework of CP violating NMSSM. The rephasing invariant 
combinations of the physical CP phases contributing to the EDMs are {(j)\ + (f)Af), {4>'x + (pi)y 
and {(p'x — (f)'i^), with (pAf and (pi denoting the CP phases of the soft trilinear parameters 
Af and the three gaugino mass parameters Mj=i^2,3, respectively. Unlike the MSSM the 
non- vanishing CP phase (0a~0k) could result in significant CP- violating mixing among the 
neutral Higgs bosons at tree level even when sin(0';!^ + (pAj) = sin(0';!^ + 0j) = 0. Throughout 
this work, we have taken a convention of (p'^ = 0. 

The MSSM CP phases (pAf and (pi are generally known to be strongly constrained 
by the non-observation of the Thallium, neutron, and Mercury EDMs, without invoking 
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cancellations among various comparable contributions, if the sfermions are within the reach 
of the LHC. On the other hand, the matter-dominated Universe requires sources of CP 
violation other than the single Kobayashi-Maskawa phase §. In this work, we concentrate 
on the CP phase 0'^, which only exist in the NMSSM framework, by taking sin(0^^) = 
sin(0j) = 0. 

One of the attractive features of the NMSSM, compared to the MSSM, might be that 
the mechanism of EWBG could be realized in a more natural setting. In analyzing the 
EDM constraint on the CP phase 0'^, we have taken a scenario in which a first-order phase 
transition is presumed to occur. The strong enough first-order phase transition is essential 
to the EWBG. 

Previously, the constraint on 0'^ from the neutron and electron (equivalently. Thal- 
lium) EDMs had been considered but only the one-loop neutralino contributions were taken 
into account. We check that our one-loop results agree well with the previous ones, showing 
no strong constraints on 0'^ especially when the sfermions of the first two generations are 
heavier than ~ 300 GeV. In this work, we have further considered the constraint from the 
non-observation of Mercury EDM and found the similar results. 

In addition to the one-loop contributions, we have taken account of the higher-order 
corrections to the EDMs and CEDMs of the light quarks and electron, the dimension- 
six Weinberg and the Higgs-exchange four-fermion operators. We note that most of these 
operators are generated due to the CP- violating neutral Higgs-boson mixing induced by the 
CP phase (f)'^. The two-loop contributions and, especially, the coefficient of the dimension- 
six Weinberg operator start to dominate when the sfermions of the first two generations 
are heavier than ~ 300 GeV. We found that they can saturate the current bound on the 
neutron EDM and can go over that on the Mercury EDM. For the Mercury EDM, we have 
included the uncertainties in the calculation of the Schiff-moment-induced term and found 
that there is still a room to have the maximal CP phase 0'^ ~ 90°. We have also shown 
that the large CP phase 0'^ ~ 90° can be easily probed in the proposed future experiments 
searching for the EDMs of deuteron and the ^^^Ra atom and it might be connected to the 
EWBG, providing a new mechanism for it [HI]. 

Furthermore, we offer a few more comments as follows before we close. 

1. The new CP phase that we considered here can be applied to a number of low-energy 
phenomenologies, such as CP asymmetries in B mesons and K mesons. 

2. In our previous work [36], we have imposed the following constraints to the parameter 
space: the LEP limits, the global minimum condition, and the positivity of the square 
of the Higgs-boson masses. In the present work, we have further constrained the 

^See Ref. [55] for the recently suggested bino-driven EWBG scenario exploiting the CP phase of the 
bino mass parameter, 0i, to account for successful EWBG without inducing large EDMs. 
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parameter space required by the EDM constraints. Therefore, by combining all these 
constraints we can explore further phenomenologies, including the Higgs physics at 
the LHC, the muon anomalous magnetic moment, and electroweak baryogenesis. 

3. One of the reasons why the 2- loop BZ diagrams can dominate is that one or more of 
the neutral Higgs bosons become very light. One can imagine that contributions to 
the muon anomalous magnetic moments can also become important in comparison 
to the one- loop result. In fact, one can show that the muon anomalous magnetic 
moment and EDM are related to the real and imaginary parts of the combination of 
couplings. We will come back to this issue in a future work. 

4. We will soon make available a computer code for calculating all the couplings and 
masses of the Higgs bosons, with the parameter space restricted by all the experi- 
mental constraints (all the above mentioned ones as well as direct search limits, muon 
g — 2, etc) as we proceed further in this framework. 
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Appendices 

Throughout Appendices, we are following the conventions and notations of CPsuperH [27] . 

A Masses and mixing matrices 

Here we set up our conventions and notations of the masses and mixing matrices of neutral 
Higgs bosons, charginos, neutralinos, third-generations sfermions. 

• Neutral Higgs bosons: 

(0°, 00, 0° , a, as)l = 0^,{H,, H,, H^, H,, H,)J , (A.l) 
where O^MjjO = diag (M^^, M^^, M|^, M^J with Mh, < Mh, < Mh, < 



25 



Charginos: We adopt the convention H^^^ = H^^^y The chargino mass matrix in 

Ma V2Mw cp 



the {W-,H-) basis 



M 



c 



(A.2) 



V2 



is diagonahzed by two different unitary matrices CrJ^cC]^ = diag{m~±, m-^:}, where 
m~± < 'm~±. The chargino mixing matrices {CL)ia and {Cjijia relate the electroweak 
eigenstates to the mass eigenstates, via 

X-R = {Cn)U-R. XaH = {W-,H-)l. (A.3) 

We use the following abbreviations throughout this paper: sp = sin/3, cp = cos (3, 
tjs = tan/3, S2/3 = sin 2/3, C2/3 = cos 2/3, sw = sin6'iy, cw = cos6'iy, etc. 

Neutralinos : The symmetric neutralino mass matrix in the {B, W^, H^, S) basis 
is given by 

MzSfjSw 
-Mzspcw 



M 



N 



(Ml -Mzcpsw 
M2 Mzcpcw 









V2 ^ 



V2 ^ 



\ 



vse 



(A.4) 



This neutralino mass matrix is diagonahzed by a unitary matrix A^: N*AiNN'' = 
diag (m~o, m~o, m~o, m~o, m~o) with m~o < m~o < m-o < m~o < m~o. The neutralino 

" ^ Al ' A2 A3 A4 A5' Al A2 A3 A4 A5 

mixing matrix relates the electroweak eigenstates to the mass eigenstates via 

{B, W', Hi Hi S)l = N:^{xI Xl Xl Xl Xl)^ ■ (A.5) 

Stops, sbottoms, staus and tau sneutrino: At the tree level, the Yukawa couplings 
are given by 



h, 



mi 



rrid 



Mi 



(A.6) 

The stop and sbottom mass matrices may conveniently be written in the {q^, qji) 
basis as 

M% + m\ + C2^M| (T*^ - Q^sl,) h\vMl^-'' ' ^e^(^-+'^) cot /3)/v^ 



%/2 



e-i(<^A+</') cot/3)/V2 



U3 



mr 



C2pMl Qts^ 



M 



V2 ^ 



Kvd{Al - ^e'('^^+^+^) tan 



tan/3)/V2 + ml 



(A. 
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with T^]^ = —T^i = 1/2, Qt = 2/3, Qb = —1/3, and hq is the Yukawa couphng of the 
quark q. On the other hand, the stau mass matrix is written in the (tl, tr) basis as 

+ml + C2pMl {si, - 1/2) KvMl - ^e^(^^+^+'^) tan/3)/v^ 
^--y h,v,{A, - ^e-(<^^+^+^) tan/3)/v^ M| + ml - c^pMl s^, 

and the mass of the tau sneutrino Ur is simply = + |c2/3M|, as it has no 

right-handed counterpart in the NMSSM as in the MSSM. 

The 2x2 sfermion mass matrix Aij- for f = t,b and r is diagonahzed by a unitary 
matrix U^: U^'^ Ai'lU^ = diag(m~,m~) with m~ < m~. The mixing matrix 

/l_ /2 /l /2 _ 

relates the electroweak eigenstates fL,R to the mass eigenstates /i,2, via 



f 



(A.9) 



B Higgs-boson couplings to sfermions 

Here we present the couplings of the neutral and charged Higgs bosons to squarks and 
sleptons in the weak basis. 



The neutral Higgs couplings to sfermions: 

— In the {bi, hRj basis, the neutral Higgs couplings to the sbottoms: 



b*b 



b*b 



L/i* A* 

\hb?vcp + \g''^vci3 



1 12 



\hl\X\vsj3 e 



b*b 



■pag b*b 



1 

71 



^ ^hl(spAl + cp\^e^^^-- 



■e+Lp) 
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— In the (tLytR) basis, the neutral Higgs couphngs to the stops 



Y<j>% t*t 



1 12 











t*t 



ilh*\X\vcpe'^^^+'^^ 
In the (jL, tr) basis, the neutral Higgs couplings to the staus: 



1 /2 

-y VSp 



T T 



V2 





/ 







r r 



2 ,.;|A|t;s;3e^(<^^+^+^) 





The neutral Higgs couplings to the sneutrinos 

r^-?- = -\ {g' + g") vcp, r^-^- = i + 

and the other couplings are vanishing. 
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The charged Higgs couphngs to sfermions: 
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Figure 8: The EDMs and CEDMs of the electron and the hght quarks together with their 
constituent contributions taking |A| = 0.81,^^1 = 0.08, I^a] = 575 GeV, and \A^\ = 110 
GeV. The other parameters are fixed as in Eq. (ISTl) . Each frame shows \df /e\ (upper left), 



10^ 



102^ 



10 



1 _ 



Ph 



1 



^ 10 



10 



10 



10 



-1 _ 



-2 



-3 _ 



-4 



Tl 
n(QCD 
n(CQM 
n(PQM 



0'. = 90 



o 




Hg(I) 
Hg(II) 
Hg(III) 
Hg(IV) 
D 
Ra 



I I I I \ \ rr^ \ \ I I I I I \ \ \ rr^ i i i i \ \ rr^ \ \ i i i i i 

100 200 300 400 500 600 700 800 900 1000 



M^, , [GeV] 



Figure 9: The observable EDMs taking |A| = 0.81, \k\ = 0.08, \Ax\ = 575 GeV, and 
U«| = 110 GeV as functions of Mr; = Mp, = Mr? = Mp; = My = Mt^ . We have 
fixed (j)'^ = 90° and the other parameters are taken as in Eq. fl5T]) . The hues are the same 
as those in Fig. [5l 
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